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ABSTRACT 

The  strain  response  of  metal  plates  during  welding  is 
discussed  and  current  state-of-the-art  efforts  to  analyze 
the  phenomena  are  reviewed.   Mechanical  and  physical  temp- 
erature dependent  properties;  transient  strain  and  temper- 
ature distribution  data  obtained  from  experiments  on 
6061  aluminum  alloy  in  the  T651  condition;  and  a  data 
reduction  computer  program  are  presented.   The  experiments 
were  designed  to  look  at  the  macroscopic  effects  of  welding 
upon  heat  treatable  age  hardened  alloys  while  at  the  same 
time  approximating  ship  structural  weldments. 

The  transient  strain  response  was  found  to  be  pre- 
dominantly longitudinal  but  transverse  strains  were  signifi- 
cant in  the  region  of  the  welding  arc.   It  was  also  found 
that  residual  strains  correlated  reasonably  between  similar 
experiments  with  a  variation  of  the  second  stress  deviator 
tensor  invariant  J2  representing  the  "state  of  stress" 
rather  than  by  individual  strain  components. 

Several  recommendations  are  made  concerning  continued 
experimental  investigation  aimed  at  further  development 
of  the  National  Aeronautics  and  Space  Adminstration  programs 
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I;   INTRODUCTION 

A.   Background: 

The  phenomenon  of  thermal  stresses  and  strains  resulting 
from  welding  processes  lead  to  difficulties  in  the  fabrication 
of  structures.   Non-uniform  temperature  distribution  leads 
to  residual  stresses  which  in  turn  provide  problems  from 
stress  corrosion,  buckling,  brittle  fracture,  fatigue  strength, 
and  distortion.   Depending  upon  the  type  of  material  being 
welded,  any  or  all  of  the  above  problems  may  be  important 
considerations  in  the  structure's  design  and  its  service  per- 
formance (1) .   All  contribute  to  the  reliability  of  the 
basic  structure  (2,  pp.  2).   Masubuchi  (2-5),  Klein  and 
Masubuchi  (6) ,  and  Klein  (7)  provide  interpretive  reports 
which  review  the  particulars  of  each  of  the  above  phenomena 
for  aluminum  and  steel  as  well  as  comprehensive  reviews  of 
the  literature.   Almost  all  of  the  previous  studies  refer 
to  residual  stresses  and  distortion  and  only  recently  has 
much  attention  been  given  to  the  actual  thermal  stresses 
during  welding.   This  is  due  to  the  complexity  of  the  problem 
characterized  by  large  temperature  changes  in  small  areas 
near  the  welding  arc  with  its  resulting  non-elastic  deforma- 
tion, temperature  dependency  of  the  material  properties  and/or 
phase  transformations,  and  complex  boundary  conditions 
resulting  from  conditions  of  geometry  and  multi-pass  welding 
(3,  pp.  2,  3)  . 
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B.   Thermal  Stresses  During  Welding/  Residual  Stresses  and 
Distortion: 

In  order  to  study  the  effects  of  residual  stresses 
and  distortion  performance  of  welded  structures,  it  is  first 
necessary  to  understand  the  relative  magnitude  and  distribution 
of  residual  stresses  and  distortion  in  weldments.   The 
mechanism  by  which  thermal  stresses  and  strains  are  developed 
in  welded  plates  is  best  described  by  Masubuchi  (8,3,  pp.  4-6) 
and  is  repeated  here: 

Figure  1  shows  schematically  changes  of  temperature 
and  stresses  during  welding.   A  bead-on-plate  weld  is  being 
made  along  the  x-axis.   The  welding  arc,  which  is  moving  at 
a  speed,  v,  is  presently  located  at  the  origin,  0,  as  shown 
in  Figure  la. 

Figure  lb  shows  temperature  distribution  along  several 
cross  section.   Along  Section  A-A,  which  is  ahead  of  the 
welding  arc,  the  temperature  change  due  to  welding,  AT,  is 
almost  zero.   Along  Section  B-B,  which  crosses  the  welding 
arc,  the  temperature  distribution  is  very  steep.   Along 
Section  C-C,  which  is  some  distance  behind  the  welding  arc, 
the  distribution  of  temperature  change  is  as  shown  in 
Figure  lb- 3.   Along  Section  D-D,  which  is  very  far  from  the 
welding  arc,  the  temperature  change  due  to  welding  again 
diminishes. 

Figure  lc  shows  the  distribution  of  stresses  along 
these  sections  in  the  x-direction,  a  .   Stress  in  the 

A 
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y-direction,  a  .  and  shearing  stress,  Txv'  also  exists 
in  a  two-dimensional  stress  field.* 

Along  Section  A-A,  thermal  stresses  due  to  welding  are 
almost  zero.   The  stress  distribution  along  Section  B-B  is 
shown  in  Figure  lc-2.   Stresses  in  regions  underneath  the 
welding  arc  are  close  to  zero,  because  molten  metal  does 
not  support  loads.   Stresses  in  regions  somewhat  away  from 
the  arc  are  compressive/  because  the  expansion  of  these 
areas  is  restrained  by  surrounding  metal  that  is  at  lower 
temperatures.   Since  the  temperatures  of  these  areas  are 
quite  high  and  the  yield  strength  of  the  material  is  low, 
stresses  in  these  areas  are  as  high  as  the  yield  strength 
of  material  at  corresponding  temperatures.   The  magnitude 
of  compressive  stress  passes  through  a  maximum  with  increas- 
ing distance  from  the  weld  or  with  decreasing  temperature. 
However,  stresses  in  areas  away  from  the  weld  are  tensile 
and  balance  with  compressive  stresses  in  areas  near  the 
weld.   In  other  words, 


fax      •   dY  =  0  (1) 

across  Section  B-B.**   Thus,  the  stress  distribution  along 

In  general  three-dimensional  stress  field,  six  stress 

components,  CT  ,  a    ,  a    ,    T   ,  T   ,  T    exist. 
ux  uy'  °z   Lxy   Lzy   lzx 
** 

Equation  (1)  neglects  the  effect  of  0      and  T   on  the 
equilibrium  condition.  y      Xy 
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Section  B-B  is  as  shown  in  Figure  lc-2. 

Stresses  are  distributed  along  Section  C-C  as  shown  in 
Figure  lc-3.   Since  the  weld  metal  and  base  metal  regions 
near  the  weld  have  cooled,  they  try  to  shrink  causing 
tensile  stresses  in  regions  close  to  the  weld.   As  the 
distance  from  the  weld  increases,  the  stresses  first  change 
to  compressive  and  then  become  tensile. 

Figure  lc-4  shows  the  stress  distribution  along  Section 
D-D.   High  tensile  stresses  are  produced  in  regions  near 
the  weld,  while  compressive  stresses  are  produced  in  regions 
away  from  the  weld.   The  distribution  of  residual  stresses 
that  remain  after  welding  is  completed  as  is  shown  in  the 
figure. 

The  cross-hatched  area,  M-M,  in  Figure  la  shows  the 
region  where  plastic  deformation  occurs  during  the  welding 
thermal  cycle.   The  ellipse  near  the  origin  0  indicates  the 
region  where  the  metal  is  melted.   The  region  outside  the 
cross-hatched  area  remains  elastic  during  the  entire  weld- 
ing thermal  cycle. 

As  shown  in  Figure  1,  thermal  stresses  during  welding 
are  produced  by  a  complex  mechanism  which  involves  plastic 
deformation  at  a  wide  range  of  temperatures  from  room  temp- 
erature up  to  the  melting  temperature.   Because  of  the 
difficulty  in  analyzing  plastic  deformation,  especially  at 
elevated  temperatures,  mathematical  analyses  were  limited 
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4. Section   D-D 


b.  Temperature  Change         c.  Stress  <r 


Figure  1     Schematic   Representation  of  Changes  of  Temperature  and  Sti 
During     Wetting     (3.8). 
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for  very  simple  cases  such  as  spot  welding. 

However,  on  the  basis  of  recent  developments  in  com- 
puter technology,  it  appears  that  a  technological  break- 
through as  far  as  the  analysis  of  thermal  stresses  during 
welding  are  concerned  is  not  too  far  away. 

C.   Previous  Investigations; 

The  closest  approach  to  a  breakthrough  to  date  has 
been  the  one-dimensional  strip  theory  approach  first  pre- 
sented by  Tall  (9)  and  extended  by  Masubuchi  and  associates 
at  Battelle  Memorial  Institute  (10)  and  Massachusetts 
Institute  of  Technology  (7,  11,  12)  under  the  support  of 
NASA's  George  C.  Marshall  Space  Flight  Center.   These 
studies  have  been  reviewed  in  detail  in  references  (3-  6) 
and  will  not  be  repeated  here. 

A  great  deal  of  effort  in  the  last  several  years  has 
been  devoted  to  developing  finite  element  programs  which 
will  extend  one-dimensional  analysis  to  two  and  three 
dimensions  (13  -  20) .   Nearly  all  have  met  some  qualitative 
success  but  in  general,  quantitative  results  have  been 
lacking.   The  finite  element  method  utilizes  an  approximate 
model  to  derive  a  set  of  equations  which  is  then  solved 
exactly.   The  flexibility  of  the  finite  element  method  will 
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allow  the  extension  of  the  flat  plate  problem  into  more 
complex  and  real  life  geometries.   It  provides  a  complete 
strain  picture,  including  longitudinal,  transverse,  and 
shear  strains.   It  has  the  capability  to  model  the  restraint 
conditions  in  the  actual  structure  and  weld.   Its  greatest 
drawback  is  high  expense  from  considerable  computer  time 
for  large  degree  of  freedom  problems  and  the  post  analysis 
of  the  vast  quantities  of  data  the  technique  generates. 
Interactive  graphics  will  greatly  aid  this  latter  problem 
in  the  future. 

Masubuchi  and  Iwaki  (13)  and  Masubuchi  and  Andrews  (14) 
coupled  the  elastic-plastic  analysis  developed  by  Wang  (15) 
with  the  thermal  loading  calculations  for  moving  heat 
sources  previously  utilized  by  the  one-dimensional  programs. 
Reference  (13)  was  qualitative  and  reference  (14)  provided 
some  quantitative  results  with  course  meshes  but  appeared  to 
have  transition  element  instability  when  a  smaller  mesh  was 
utilized.   In  this  situation,  an  element  which  reaches 
yield  stress  is  considered  a  transition  element.   However, 
in  the  next  iteration,  the  equivalent  stress  from  that 
element  drops  below  yield  and  the  element  becomes  elastic 
again.   The  whole  situation  jumps  between  the  two  sets  of 
values.   In  order  to  overcome  this  difficulty,  a  scheme 
was  devised  to  quarter  the  time  step   in  which  this  insta- 
bility occurred.   This  worked  for  a  large  number  of  cases 
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but  not  all  cases.   Sometimes  increasing  the  time  step 
also  eliminated  the  instability.   Reference  (10)  may  offer 
a  way  to  correct  this  problem.   Yamada  (17)  seems  to  have 
met  with  more  success  dealing  with  thermal  stresses  due  to 
rapid  heating  by  controlling  the  number  of  elements  which 
yield. 

Hibbit  and  Marcal  (18)  offer  the  most  analytic  of  the 
approaches  to  date.   Their  model  treats   the  weld  process 
as  a  thermo-mechanical  problem.   A  finite  element  formula- 
tion derived  from  the  uncoupled  thermal  and  mechanical 
energy  balances  forms  the  basis  of  the  model.   It  attempts 
to  deal  with  material  non-linearity  due  to  temperature 
dependence  of  thermal  properties  and  in  the  fusion  problem 
where  material  phase  change  is  accompanied  by  a  latent  heat 
effect.   It  includes  radiation  as  a  cooldown  mechanism  and 
finite  strain  effects  during  isothermal  loading.   However, 
there  was  little  agreement  with  experimentally  measured 
residual  stresses  and  they  concluded  that  their  finite 
element  model  did  not  include  significant  material  behavior. 

Shinn  (19)  approaches  the  two-dimensional  distortion 
of  a  panel  structure  due  to  welding  with  the  assumption 
of  elastic  deformation  during  the  welding  process.   Com- 
putations were  carried  out  using  the  one-dimensional 
experimental  values  of  unconstrained  angular  changes  along 
the  welded  edge  and  its  equivalent  constrained  welding 
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moment  as  an  input  to  the  computer  program.   The  results 
were  not  completely  successful  because  of  the  elastic 
assumption  and  the  questionable  accuracy  of  the  experimental 
data  input  into  the  program.   It  is  of  the  opinion  of  this 
author  that  this  technique,  on  the  basis  of  current  tech- 
nology/ will  provide  the  most  fruitful  results  in  the 
immediate  future  in  terms  of  financial  and  manpower  invest- 
ments.  The  basic  welding  problem  is  highly  non-linear 
and  may  be  too  complex  to  solve  in  a  completely  analytical 
form.   A  purely  empirical  approach  does  not  appear  very 
fruitful,  either.   It  is  not  a  simple  task  to  determine 
thermal  stresses  in  small  regions  heated  to  high  temperatures. 
Without  proper  analysis,  it  would  be  difficult  to  adequately 
interpret  experimental  data. 

The  empirical-analytic  technique,  however,  is  only  as 
good  as  experimental  data  provided.   At  present,  material 
and  physical  properties  at  the  range  of  temperature  from 
room  to  melting,  even  for  the  most  common  of  structural 
materials,  are  sorely  lacking  in  the  literature.   More 
effort  must  be  expended  in  providing  this  critical  information, 
The  approximations  utilized  in  the  interim  will  be  the 
controlling  feature  of  this  computational  technique. 
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D.   Aim  and  Purpose  of  the  Present  Study: 

Muraki  and  Masubuchi  are  presently  developing  a  new 
elasto-plastic  finite  element  computer  program  at  Massachusetts 
Institute  of  Technology  along  the  lines  described  by  Yamada 
(17)  with  temperature  distribution  adopted  from  reference  (20) . 
It  is  the  purpose  of  this  report  to  provide  experimental 
data  on  a  heat  treatable  tempered  aluminum  alloy  to  verify 
this  program.   Electric  resistance  thermocouples/  three 
element  strain  gage  rosetts,  single  element  strain  gages, 
and  an  extensioneter  will  be  utilized  to  obtain  transient 
and  residual  strain  and  displacement  data  in  the  center  and 
near  the  edge  of  30x18  inch  panels  for  an  automatic  GMA 
welding  process.   Experiments  will  be  performed  for  bead-on- 
plate  and  butt  welding  processes  under  similar  heat  inputs 
to  verify: 

1.  that  experiments  may  be  duplicated  under  reasonably 
similar  conditions  to  give  repeatable  results; 

2.  to  compare  butt  and  bead-on-plate  welds  for  similar- 
ities and  differences; 

3.  to  compare  strains  in  both  top  and  bottom  of  the 
plates  to  study  the  importance  of  bending  strains  as  com- 
pared with  transverse  and  longitudinal  strains; 

4.  to  observe  transient  principal  strain  magnitudes 
and  directions; 

5.  to  measure  residual  stress/strain  distributions 
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in  the  longitudinal  and  transverse  directions; 

6.  to  measure  the  deflection  of  the  two  half  plates 
in  a  butt  welding  process  under  tack  welded  conditions; 

7.  to  present  longitudinal  and  transverse  strains  for 
computational  technique  verification,  and; 

8.  to  collect  and  present  in  usable  form  available 
physical  and  mechanical  temperature  dependent  data  for 
use  in  the  computer  analysis. 

Standard  strain  gages  will  be  utilized  which  may  not 
exceed  400°F.   Since  the  material  properties  in  the  heat 
affected  zone  (HAZ)  are  not  fully  known  or  understood  and 
the  high  temperature  strain  gages  have  proven  unreliable 
(6,  7)  in  part  because  of  temperature  compensation  data  and 
in  part  by  the  lack  of  material  properties  data  within  the 
HAZ,  strain  measurements  will  be  beyond  the  HAZ. 
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II;   MATERIAL  BEHAVIOR 

A.  General; 

A  new  finite  element  program  to  study  residual  stresses 
and  metal  movement  is  presently  being  developed  by  M.I.T.'s 
Department  of  Ocean  Engineering  concurrent  with  this  thesis. 
Since  an  attempt  is  being  made  here  to  provide  material 
properties  for  and  consistent  with  this  new  program,  a 
brief  and  simple  formulation  of  the  tangent  stiffness 
technique,  incremental  stress  strain  relationship,  and 
thermal  loading  terms  is  provided  for  the  sole  purpose  of 
demonstrating  the  importance  of  temperature  dependency.   Then, 
material  properties  consistent  with  this  formulation  are 
presented  for  6061  aluminum  alloy  in  the  T6  and  T651  condi- 
tion, the  material  utilized  for  the  experiments  presented 
in  Chapter  III. 

B.  Mathematical  Formulation; 


The  most  common  finite  element  approach  to  plasticity 
problems  in  welding  has  been  the  Tangent  Stiffness  Approach. 
The  problem  is  first  solved  elastically  for  the  general 
loads.   Then  the  elastic  limit  loads  are  determined  by 
scaling  down  the  given  loads  until  the  element  with  the 
maximum  equivalent  stress  is  just  reached.   The  differences 
between  the  given  loads  and  the  elastic  limit  loads  are  then 
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Figure  2:  Equivalent  Stress  and  Strain  (14, pp  20) 
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Figure  3:   Iteration  Process  to  Bring  Back  the  Stress  Point 
on  to  the  <!f- €.  curve  with  H*.  The  element  in  the  tran- 
sition region  starts  out  at  point  B  as  elastic.  The  elastic 
increment  BC  of  the  element  for  the  prescribed  load  increment 
dP  is  first  calculated.  The  stress  point  corresponding  to 
C  should  be  D  on  the  stress-strain  curve.  Using  the  modulus 
corresponding  to  BD,  the  increment  BF  for  dp  is  next  calcu- 
lated. This  process  may  be  repeated  until  the  solution 
converges   (17,  pp  300,303). 
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divided  into  a  predetermined  number  of  loading  increments. 
The  problem  is  then  solved  including  plasticity  for  each 
element.   When  an  element  becomes  plastic,  the  tangent 
stiffness  approach  is  used  to  determine  an  effective  elasto- 
plastic  stiffness  for  that  element.   In  this  approach, 
the  Prandtl-Reuss  relation  and  the  plastic  modulus  are 
combined  with  the  elastic  stress-strain  relations  to  produce 
an  incremental  stress-strain  relation  which  includes  the 
effects  of  plasticity.   These  incremental  relations  depend 
upon  the  instantaneous  stress  levels  in  each  element  (14, 
pp.  4  -  5) . 

In  general, 

da   =   E.  de  (2) 


where  E   is  the  tangent  modulus  coefficient.   E.  has  an 
elastic  term  and  a  plastic  term.   The  plastic  term  is  zero 
if  no  yielding  takes  place.   Then,  E,  may  be  used  to  estab- 


lish the  tangent  stiffness  array  K   for  each  element, 

<-3 


K 
n 


Bt  E,_   B  d  (Volume)  (3) 

~n  ^t  n 


Volume 

where  B  defines  the  strain  displacement  relationships. 
Thus,  where  loading  is  applied  in  increments,  the  structure 
stiffness  array  may  be  found  at  each  stage  of  loading  by 
evaluating  each  element's  tangent  stiffness  (21,  pp.  12  -  13) 
The  incremental  pseudo-thermal  loading  term  likewise 
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becomes 


bJ;  E.  de.  d  (Volume)  (4) 

n  y-s3*    -  -^  - 

Volume*"* 


where  de   is  the  thermal  strain  caused  by  the  welding 
process.   Summed  over  the  entire  structure,  one  obtains: 


K  dU   =   dP  (5) 


dU  =   K_1  dP  (6) 


de   =   bF  dU  (7) 


da   =  Et  Bt   dU  (8) 


To  determine  the  tangent  modulus  coefficient  related 
to  equation  (2),  one  must  consider  the  following  relations:* 

Prandtl-Reuss  Flow  Rule: 

deP   =  |£  diP  (9) 

Yielding  Surface  Behavior: 

-  T 
da  =  |2.    da  (10) 

aO  ^v^ 

Material  Behavior: 


a     =  a0  +  eP  F   =  a0  +  H'  ep        (11) 

da  =  e£  de"   =  h'  de„  (12) 

t   p  P 

Utilizing  equations  (9)  -  (12)  and  recalling  that 


The  p  and  e  superscripts  refer  to  plastic  and  elastic 
respectively. 
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da   =   Dde   =   D  (de  -  de*)  (13) 

where  de  is  the  total  incremental  strain  and  D  is  given  in 

Equation  (24)  for  the  plane  stress  situation,  one  arrives 

at  the  following  relationship: 


"  D  H 


8a 


H1  +   l£  T  D  |£ 
9a   ~  3i 


de 


(14) 


This  reduces  to  equation  (2)  where  the  quantity  in  j  V 
is  the  tangent  modulus  coefficient.  This  is  consistent 
with  reference  (17,  pp.  310) . 

If  we  assume  plane  stress  conditions  and  consider  the 
equivalent  stress  for  Von  Mises1  isotropic  material, 


—  2       2  2         2 

a   =  a    -     a  o        +     a    +3t 

x       x  y      y        xy 


(15) 


Differentiating, 


da 


=  —   (2a   -  a 
2a  L    x 


y)  dax  +   (2ay 


This  may  be  rewritten  in  matrix  form  as 


*-     1   „T 
da  =  —  a 


l-J  0 

"I   *   ° 

0   0   3 


da 


or 


-     IT 

da  =  -  a  F  da 


a  )  da   +  6t  dx 
x7    y       xy  xy 

(16) 


(17) 


(18) 
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where 


xy 


(19) 


da  = 


J 


da, 

d0; 

df 


xy 


(20) 


Comparing  equations  (10)  and  (18) ,  we  see  that 

—  T 

8(7      1   T  _ 
__   =   _  0      p 

do  ~    ~ 

s+*  a 


(21) 


Transposing, 


8a  1  PT„  1  „  n 
-~ —  =  —  Fa  =  —  Fa 
3a     -  ~^  v~*~ 


(22) 


Substituting  the  previous  relationships  into  equation  (14) 
we  arrive  at 


da 


D  F  a   a   F  D 

—  2        T 

~    a  H'  +  a  F  D  F  a 


de 


(23) 


for  a  two  dimensional,  plane  stress,  isotropic  case,  with 
the  Von  Mises'  yield  criterion  where 


D  = 


1  -  v' 


1   v 
0    1 


0 

0 

1-v 


F   = 


1  -|   0 


-i   i    ° 


-2        2 

a        =      a 
x 


(24) 


(25) 


-     o  o       +  a„ 
x  y      y 


+   3t 


xy 


(15) 


28 


Equation  (23)  is  but  a  simple  example  of  the  more 
complex  general  case.   The  important  point  to  note  here  is 
that  E,  o,  and  H1  are  functions  of  temperature  and  this 
makes  quite  complex  even  the  most  simple  of  welding  problem 
examples. 

In  the  previous  study  at  the  Massachusetts  Institute 
of  Technology,  Andrews  (14,  pp.  25)  assumed  the  simplest 
model  possible  for  the  nodal  force  matrix,  i.e.,  that  of 
a  constant  temperature  of  each  element.   In  this  case,  the 
nodal  force  matrix  for  each  increment  becomes 

dP„   =       (BT„  E.  de.  )  d  (Volume)        (4) 


Volume 


det   =   AT 


where 


AT  =   average  temperature  change  of  element 
a  =   coefficient  of  linear  thermal  expansion 

If  the  element  in  question  is  plastic,  the  tangent  stiffness 
matrix  is  used.   Masubuchi  and  Iwaki  (13,  pp.  9)  assumed 
a  much  more  complex  loading  function  which  accounted  for 
the  equivalent  nodal  forces  due  to  the  dependency  of  the 
coefficient  of  thermal  expansion  and  Young's  modulus  on 
temperature.   The  plastic  terms  additionally  took  the 
temperature  dependency  of  the  yield  surface  into  account. 
Yielding  is  complex.   It  can  occur  during  the  heating  stage 
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and  also  in  the  subsequent  cooling  stage  since  some 
elements  can  be  stressed  to  their  yield  point  in  the 
reverse  direction.   Yamada  (17,  pp.  298)  assumed  the  same 
simplified  type  of  expression  as  equation  (4)  in  his 
study  of  thermal  stresses  due  to  rapid  heating.   As 
in  the  case  of  the  incremental  stress-strain  relation,  the 
important  point  to  note  is  that  in  varying  degrees  of 
difficulty,  the  loading  term  is  dependent  upon  E,  a,  and 
a  which  are  in  turn  dependent  upon  temperature. 

The  temperature  distribution,  uncoupled  from  the 
mechanical  problem  just  described,  may  be  approximated  by 
classical,  finite  difference,  or  finite  element  techniques 
(9-11,  13,  14,  18,  20,  24-30).   In  terms  of  the  conductiv- 
ity tensor,  k,  the  general  representation  of  the  aniso- 
tropic heterogeneous  continuum  is 

pc  ||  =   V  •  (k  •  VT)   +   Um  (25) 

where  p  is  the  density,  c  is  the  specific  heat,  k  is  the 
thermal  conductivity,  and  U,M  is  the  heat  generation  per 
unit  volume  (23,  pp.  29,  44-45).   Fourier's  law  applied 
to  the  boundary  conditions  provides 

q  =   -k  •  VT  (26) 
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where  q  is  the  heat  flux  normal  to  the  boundary  surface. 
If  the  continuum  is  isotropic,  these  vector  equations 
reduce  to  the  familiar 


dT 
pc  gi  =   V  •  (kVT)   +   Um  (27) 


q   =   -kVT  (28) 

Equation  (27)  may  be  reduced  to 


dT  2 

pc  g|  =   Vk  •  VT   +   kV  T   +   U"'  (29) 


If  k  is  a  function  of  space  only,  equation  (29)  is 
linear.   On  the  other  hand,  where  <  depends  on  temperature 
alone, 
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equation  (29)  may  be  modified  to: 

PC  ||  =   ||  (VT)2   +   kV2T  +   u'"  (30) 


which  is  nonlinear.   For  homogeneous  isotropic  continua  k 
is  constant  and  equation  (29)  reduces  to 

dT       2      Um 

dt  pc  (31) 

where  a  is  the  thermal  diffusivity.   Since  p,  c,  and  k 
are  functions  of  temperature,  the  welding  problem  becomes 
highly  nonlinear.   Coupled  with  phase  changes,  this 
problem  becomes  even  more  complex. 

C.   Material  Properties: 

As  discussed  in  Section  B,  tabulation  of  the  material 
properties  as  a  function  of  temperature  is  of  prime  impor- 
tance.  This  information  is  difficult  to  obtain  and  often 
simply  is  not  available.   In  this  experimental  work,  alloy 
6061  in  the  T651  Condition  is  utilized.   It  is  one  of  the 
most  versatile  of  the  wrought,  heat  treatable  aluminum 
alloys.   An  understanding  of  the  behavior  of  this  common, 
relatively  inexpensive  material  during  welding  and  heat 
treating  may  also  contribute  to  the  understanding  of  the 
behavior  of  quenched  and  tempered  steels  such  as  HY-80  and 
HY-130  currently  in  use  today.   6061  is  readily  welded  by 
all  methods  and  has  excellent  weldability  characteristics. 
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The  plate  used  in  this  experiment  has  a  Federal  Specifi- 
cation of  QQA-2 50/11. 
1.   Composition  (31,  Code  3206,  pp.  1) 


Percent 

Min.      Max. 

Chromium 

0.15 

0.35 

Copper 

0.15 

0-40 

Iron 

- 

0.7 

Magnesium 

0.8 

1.2 

Manganese 

- 

0.15 

Silicon 

0.40 

0.8 

Titanium 

- 

0.15 

Zinc 

- 

0.25 

Other  impurities 

each 

- 

0.05 

total 

- 

0.15 

Aluminum 

Balance 

TABLE  1:   Composition  of  6061  Aluminum  Alloy 

2.   Heat  Treatment  (31,  Code  3206,  pp.  1) 

Annealed  (0  Condition) :   Heat  to  775°F  for  2  to  6 
hours,  cool  at  50 °F  per  hour  maximum  to  500 °F.   The  rate 
of  subsequent  cooling  is  unimportant. 

T4  Condition:  Solution  heat  treat  to  970°F,  water 
quench,  and  naturally  age  at  room  temperature  to  a  sub- 
stantially stable  condition. 

T6  Condition;   Solution  heat  treat  to  970 °F,  water 
quench,  artificially  age  by  precipitation  heat  treatment 
at  320°F  for  16  to  20  hours  or  350°F  6  to  10  hours  and 
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air  cool. 

T651/T451  Condition:   Same  as  the  T6/T4  Condition  followed 
by  a  1  1/2  to  3  percent  permanent  set  stretch  stress  relief 
prior  to  any  precipitation  heat  treatment. 

The  soft,  as  quenched  0  Condition  can  be  preserved 
by  refrigeration  in  order  to  minimize  springback  and 
increase  the  general  ease  of  forming  operations.   Maximum 
holding  times  which  will  preserve  the  formability  of  this 
condition  are  shown  in  Table  2  (31,  Code  3206,  pp.  3). 


Temperature  °F 

R.T. 

32 

20 

Time 

2  hours 

2  days 

7  days  + 

TABLE  2:   Holding  Times  for  6061  in  0  Condition 


3.   Specified  Mechanical  Properties  (31,  Code  3206,  pp.  2) 
TABLE  3:   Specified  Mechanical  Properties. 


Condition 

0 

T4,  T451 

T6,  T651 

Thickness  -  in. 

0.250 

0.250 

0.250 

Ultimate  Stress  in 
Tension  min.  ksi 

22 

30 

42 

0.2%  Offset  Yield 
Stress  in  Tension 
min  -  ksi 

12 

16 

35 

Ultimate  Elongation 
2  inch  min  -  % 

18 

16  -  18 

6-10 

Ultimate  Shear  Stress 
Typical  -  ksi 

12 

24 

30 

34 


4.   Mechanical  Properties : 

In  welding  aluminum,   it  must  be  remembered  that 
the  metal  in  the  area  of  the  weld  will  be  in  an  annealed 
condition  after  welding.   There  will  be  a  corresponding 
loss  of  strength  unless  heat  treatment  is  done  after 
welding.   Welding  does  not  always  reduce  the  strength 
of  the  heat  treated  alloys  to  that  of  the  fully  annealed 
condition  because  there  is  some  air  quenching  of  the  metal 
as  it  cools  (32,  pp.  6.33).   Table  4  (32,  pp.  6.34) 
demonstrates  this  effect  for  the  type  of  filler  metal 
utilized  in  this  experiment. 


Filler 
Wire 

Specified  minimum 
tensile  strength  of 
base  plate  in  ksi. 

Average  tensile 
strength  across 
weld  in  ksi. 

Average  bend 
free  elonga- 
tion % 

Not  heat  treated  after  welding 

4043 

42.0 

27.2 

16.0 

Heat  treated  and  aged  after  welding 

4043 

42.0 

43.5 

11.0 

TABLE  4:   Strength  and  Ductility  of  Welded  Butt  Joints  in 
Aluminum  (TIG  and  GMA  with  argon) . 


Brungraber  and  Nelson  (33)  report  that  6061  alloy  may 
be  heated  up  to  550° F  by  butt  welding  without  appreciably 
affecting  the  mechanical  properties  of  the  base  plate. 
They  further  report  that  for  thicknesses  less  than  one 
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inch,  the  distance  from  the  center  of  the  weld  to  the 
edge  of  the  heat  affected  zone  (HAZ)  may  be  predicted 
as : 

bh  =   0.6t   +   0.00059  H  (32) 

where  E  has  units  of  volts,  I  amperes,  V  inches  per 
minute,  and  t  (thickness)  inches.   As  an  example,  for 

E  =  20.,  I  =  240.,  t  =  .25  and  V  =  32.16,  El/Vt  =  597 

2 
watt  -  min/in  ,  and  b,  =  .50  inch.   The  extent  of  the 

fusion  zone  may  depend  somewhat  on  the  parameter  El/Vt 

but  it  is  governed  primarily  by  the  thickness  t  of  the 

plate  and  the  geometry  of  the  edge  preparation. 

Figures  4  through  10  provide  mechanical  properties 

at  both  room  temperature  and  elevated  temperatures.   Figure 

4  provides  typical  tension  and  compression  stress  strain 

curves  for  the  experimental  plate.   Figure  5  provides 

the  average  stress  strain  curves  for  varying  temperatures 

from  room  temperature  to  700°F  based  on  a  half  hour 

exposure  to  the  temperature  in  question.   This  may  or  may 

not  be  a  valid  approximation  for  the  welding  problem,  but 

it  is  the  only  information  currently  available.   Figure  6 

shows  the  effect  of  this  temperature  on  the  0.2%  offset 

yield  stress.   This  may  be  approximated  as  the  cubic 

polynomial  given  in  Table  5.   Figure  7  demonstrates  the 

effect  of  the  half  hour  soak  temperature  on  Young's 
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modulus.   This  curve  was  generated  from  Figure  5.   Its 
cubic  approximation  is  listed  in  Table  5. 

Figure  8  shows  the  room  temperature  tangent  modulus 
at  various  stress  levels.   Its  straight  line  approximation 
is 

H'  x  10~3(ksi)   =   0.917a  (ksi)   +   38.064 

(33) 

The  highest  stress  level  shown  on  this  figure  is  40.7  ksi. 

+3 
This  equates  to  H'  =   0.746  x  10   ksi  and  may  be  utilized 

effectively  for  constant  strain  hardening.   Figure  9 

represents  an  estimate  of  the  temperature  effects  on  the 

tangent  modulus.   If  one  assumes  that 

H»  (T°F)   oc   E(T°F)  (34) 

then 

ttI    frpOrM  H'    (R.T.  )      X     E  (T°F)  Cit:\ 

H  (T  F)   -        E(R.T.)  (35) 

Masubuchi  and  Andrews  (11)  define  the  constant  strain 
hardening  parameter  for  their  one-dimensional  computer 
program  as 

H'   =  mE  (36) 

From  the  above  equations, 

m   =   H'/E   =   H' (T°F)/E(T°F)   =   0.0742        (37) 
This  is  the  first  such  value  ever  determined  for  their 
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f.p.r.4:     ««»  T«^««ti«  T«.lW^  C-»«..i»  $tr...  Stwlr 
C«r*«*  for  «<*1  1*  tl»  T*  Conditio*  (31). 
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;        .      StMtn  (Inchaa/Iwch) 
flgon  5:  SU»«  Strain  Curvoa  for  6061  In  T6  Condition  aftor 
0.5  Hour  at  Tonporatnro  (31, 32) 
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program. 

Figure  10  also  demonstrates  a  computer  algorithm 
for  the  tensile  properties  at  various  temperatures  based 
on  the  half  hour  soak.   To  determine  this  approximate 
value,  proceed  as  follows: 

a.  Determine  a  (T°F)  from  Table  5  or  Figure  6. 

b.  Determine  E(T°F)  from  Table  5  or  Figure  7. 

c.  Determine  H' (T°F)  from  equation  (35)  and  Table  5 
or  Figure  9. 

d.  Plot  o    (T°F)  at  0.2%  offset. 

e.  Draw  line  through  point  d  above  with  slope 
H'  (T°F) . 

f.  Connect  line  in  e  above  with  line  through  origin 
with  slope  E(T°F). 

5.   Physical  Properties 

Figures  11  through  14  show  the  effect  of  temperature 
on  the  coefficient  and  mean  coefficient  of  linear  expansion 
a,  density  p,  thermal  conductivity  k,  and  specific  heat 
C.   The  coefficient  of  linear  expansion  shown  in  Figure  11 
is  based  on  empirical  equations  for  high  purity  aluminum 
and  alloy  constants  provided  by  reference  (35) .   The  alloy 
constant  for  6061  TO  is  0.990,  i.e.,  0.990  a  (pure 
aluminum)  =  a  (alloy) .   The  value  of  this  constant  for 
heat  treated  tempers  is  approximately  0.015  greater. 
With  these  tempers,  application  is  restricted  to  temper- 
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atures  which  do  not  appreciably  exceed  those  used  in  the 
final  aging  treatments  (320-350°F) .   The  basic  alloy 
constant  is  limited  to  temperatures  below  600°F.   Figure  11 
is  extrapolated  past  600°F  to  the  melting  point.   Figure  12 
provides  an  estimate  of  density  provided  by  reference 
(36)  extrapolated  to  the  melting  point.   Figure  13  provides 
an  estimate  of  thermal  conductivity.   Reference  (37, 
Volume  II,  Part  2,  pp.  769)  provides  thermal  conductivity 
information  of  6000  series  aluminum  alloys  from  300  to 
650°K.   A  linear  relationship  was  plotted  parallel  to 
this  data  passing  through  the  one  data  point  provided  by 
reference  (38).   Reference  (37)  (Volume  I,  pp.  11)  listed 
the  specific  heat-temperature  relationship  for  pure  aluminum. 
A  parallel  relationship  (39,  pp.  180)  passing  through  the 
data  point  provided  by  reference  (31,  Code  3206,  pp.  1) 
resulted  in  the  estimated  data  shown  in  Figure  14. 

These  physical  values  may  be  approximated  by  the 
cubic  polynomials  given  in  Table  5  for  computer  simulation. 
Appendix  A  lists  the  computer  program  for  determining 
the  coefficients  given  in  Table  5  and  for  testing  these 
coefficients  in  10°F  increments. 

Table  6  summarizes  the  physical  and  mechanical  pro- 
perties of  6061  T6  in  the  form  utilized  by  Masubuchi  and 
Andrews  (11)  in  their  one-dimensional  program. 
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figure  II:     CMfflciMK  «f  Tlwiwl  taction  for  6061  (3$). 
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III;   PROCEDURES 

A.  Scope  of  the  Research; 

A  series  of  six  experiments  measuring  temperature 
and  strain  changes  during  welding  was  performed.   System 
models  representing  constrained  bead-on-plate  and  butt 
joints  were  constructed  from  6061  aluminum  alloy  in  the 
T6  condition.   Welding  procedures  were  utilized  which 
offered  full  penetration  with  minimum  heat  input.   This 
allowed  strain  gage  locations  as  close   to  the  weld 
line  as  possible  without  exceeding  the  4  00°F  maximum 
allowable  gage  temperature  limitation. 

B.  Selection  of  Parameters; 

The  aluminum  selected  for  this  investigation  was 
picked  because  of  its  easy  availability  and  wide  use  in 
the  marine  and  aeronautical  industries.   It  develops  its 
strength  from  its  heat  treatment.   The  0.250  inch  plate 
thickness  used  in  this  study  was  determined  by  availability, 
past  experimental  work,  ease  of  handling,  and  an  interest 
in  approaching  a  two-dimensional  problem,  yet  minimizing 
the  extent  of  bending  and  buckling.   The  joint  design 
for  the  butt  welds  used  in  this  experiment  (straight) 
was  selected  for  simplicity  and  ease  of  preparation. 
Steel  backing  plates  were  utilized  below  the  weld  and  at 
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the  clamped  ends  of  the  plates  to  provide  support  and 
distributed  clamping  and  yet  minimize  heat  conduction 
effects.  Each  specimen  was  clamped  at  its  edges  to  a 
0.50  0  inch  mild  steel  bed  plate  which  resisted  deformation. 

The  size  of  the  test  plates  was  set  at  a  nominal  18 
by  30  inches  to  provide  essentially  steady  state  conditions 
at  strain  gage  locations,  and  to  approach  the  temperature 
distribution  of  an  infinite  plate. 

The  weld  process  used  was  semi-automatic  Gas  Metal-Arc 
(GMA  or  MIG) .   It  allows  control  of  weld  variables,  reduces 
operator  error,  and  fosters  repeatability.   No  preheating 
was  necessary. 

C.   Strain  Measurement  by  Electric  Resistance  Strain  Gages: 

The  fundamental  concept  of  strain  gage  operation 
is  that  certain  conductors  exhibit  a  change  in  electrical 
resistance  with  a  change  in  strain.   Gages  designed 
according  to  this  principle  are  attached  to  test  materials 
whose  strains  are  then  monitored  by  measuring  resistance 
variations  across  the  gage.   In  the  case  of  welding 
thermal  strains,  the  observed  resistance  change,  AR,  is 
made  up  of : 

AR  =  AR..  (e  )   +   AR.(e  )   +   AR  (aT)   +   AR.  (T) 
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where 


AR,  (e  )   =   the  resistance  change  corresponding  to 

elastic  mechanical  strain,  e    ,    from  which 
stresses  can  be  computed. 

AR_ (e  )   =   the  resistance  change  corresponding  to 

plastic  mechanical  strain,  e  ,  if  it  exists 

p 

Ar  (aT)   =   the  resistance  change  corresponding  to 
temperature  induced  thermal  strain,  aT. 

AR. (T)    =   the  resistance  change  caused  by  thermo- 
electric effects  in  the  gage  itself. 

While  it  is  not  presently  possible  to  discriminate 

between  the  two  mechanical  strains,  £   and  e  ,  AR_ (aT) 

e      Pi 

and  AR. (T)  can  be  separated  out  by  physical  temperature 
compensation  of  the  strain  gage  material  coupled  with 
empirical  calculation.   For  this  investigation,  a   for  the 
gage  was  approximately  13.0  microstrain/°F.   The  AR . (T) 
and  the  difference  between  gage  a   and  base  plate  aT 
will  then  generate  an  error.   To  compute  this  error, 
a  test  gage  of  the  type  and  lot  used  in  the  experiment 
is  mounted  on  a  small  sample  of  2  024  aluminum  alloy  by 
the  gage  manufacturer.   The  sample  is  then  heated  at 
equilibrium  until  a  curve  of  "apparent  strain"  vs.  temper- 
ature is  obtained  to  the  operating  temperature  range. 
The  gage  readings  recorded  in  the  weld  experiments  can 
then  be  corrected  by  subtracting  out  the  apparent  strain 
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value  corresponding  to  temperatures  observed  at  the  gage 

location.   Figures  16  and  17  show  this  factory  provided 
information  for  the  individual  gages  and  rosettes 
respectively.   In  each  case,  the  curves  have  been 
reduced  by  the  factory  to  fourth  order  polynomials  yield- 
ing microstrain  for  ease  in  computer  calculation. 
Gage: 
EAP  =  -69.36  +  2.07T  -  2.02  x  10_2T2  +  5.39  x  10"5T3 

-  3.27  x  10"8T4  (38) 

Rosettes : 

EAP  =  -60.68  +  2.22T  -  2.40  x  10_2T2  +  7.05  x  10_5T3 

-  5.59  x  10-V  Since39) 
6061  T6  has  a  greater  coefficient  of  thermal  expansion 
than  the  2024  T4  test  plate,  a  correction  to  these  poly- 
nomials has  been  made.   This  correction  was  calculated 
as  follows: 

In  general,  T 

Thermal  Strain  = 

1 

From  reference  (35)  , 

a6061  =   C6061a  Pure  A1'   =   -99  a  pure  Al.   (41) 
a2024  =   c2024a  pure  A1'   =   -97  a  pure  Al. 


2 

a(T)  dT  (40) 

T. 


(42) 


a 

Tvi  6061  (Ar>\ 

a   pure  Al .   =   ^4JJ 

.99 


AEAP   = 


T2 

CC6061-  C2024>  "pure  Al.  (T)  dT 
Tl 


(44) 


-?.  0  +  -2 

...     =    rvVJEv'T     SHAWE    I  >J    -V»'l-4    TE:4Pfi:tArJ.W    liAo:-:    FACT'J/i 
C    O    1   ?    J   T    K     {         P    ?   'i    C    E    3    3    £    I)         n    a    r    A 
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Flgure   17:     Apparent  Strain  Correction  for  FAER-18RB-12S13-ET 
Strain  Gaga  Roaattas  Utilised   In  thla  Study. 
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AEAP  = 


T2 

(^||)  («6061CT>)  dT  (45) 

Tl 


Integrating  the  polynomial  given  in  Table  5  for  a,,.,,...  (T)  , 
0  <  T  <  400°?,  between  77°F  and  T°F  yields  the  following 
apparent  strain  correction  for  both  gages  and  rosettes: 

AEAP  =  -19.4979  +  .2463T  +  .9561  x  10"4T2  -  .7879  x  10~7T3 
+  .5303  x  10~10T4  (46) 

Equation  (46)  must  be  added  to  equations  (3  8)  and  (39) 
to  obtain  true  apparent  strain  corrections.   The  true 
apparent  strain  is  then  subtracted  from  the  actual  strain 
to  provide  true  mechanical  strains. 

Strain  gages  manufactured  in  the  form  of  rosettes 
provide  as  many  as  three  independent  readings  in  three 
directions  at  a  single  location.   This  is  enough  to 
completely  describe  the  two-dimensional  strain  state 
at  that  location  (41,  42,  43)  and  is  discussed  in  the 
next  chapter. 

D.   Description  of  Apparatus: 

1.   Test  Plates:   A  description  of  the  test 
specimens  is  summarized  in  Figures  18  through  22  which 
provide  dimensions  and  gage  locations.   Devices Q)  through 
(f2)were  dynamically  measured  during  welding.   Gages 
(^through  (%5) were  statically  measured  before  and  after 
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Figure  18:   Test  No.  1  and  Test  No.  2   (Butt  and  Bead-on-Plate 
Welds  Respectively)  Measuring  Temperature  Distribution. 


Figure  19:  Test  No.  3  (Bead-on-Plate  Weld)  Measuring  Temperature 
Distribution  and  Strains. 
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Figure  20:     Test  No.   4   (Bead-on-Plate  Weld)   Measuring  Temperature 
Distribution  and  Strains. 


&d&&  Ht-VV 


o 


-;PV 
"c  j      -: 


Figure  21:   Test  No.  5  (Butt  Weld)  Measuring  Temperature 
Distribution,  Strains,  and  Extension  of  Root  Gap. 
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Figure  22:  Test  No.  6  (Butt  Weld)  Measuring  Temperature 
Distribution,  Strains,  and  Extension  of  Root  Gap* 
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welding.   The  arrangement  of  constraining  clamps  is 
shown  in  Figure  23. 

2.   Sensors  and  Instrumentation ; 

a.  Strain  Gages.   Two  types  of  gages  were  used 

in  this  investigation:   SR-4  foil  45°  rosettes  and 
SR-4  foil  single  element  gages.   Gage  properties 
were  as  follows: 
Gage  SR-4  Rosette        SR-4  Gage 

Designation  FAER-18RB-12S13-ET    FAE-25-12S13L 

Manufacturer  BLH  Electronics    BLH  Electronics 

Grid  Length  (in.)  3/16  1/4 

Grid  Width  (in.)  0.90  0.13 

Overall  Length  (in.)  0.280  0.35 

Overall  Width  (in.)  0.540  0.13 

Temperature  Range         -50  to  +400°F     -50  to  +400°F 
Resistance  120.  120. 

Gage  Factor  2.03  2.07 

Cement  EPY-600  EPY-60  0 

Protective  Covering   BLH  Barrier-C        BLH  Barrier-C 

b.  Temperature  Sensors:   All  temperature  sensors 
used  were  Chrome 1/Alumel  thermocouples  made  from 
Leads  and  Northrup  No.  28  wires.   Each  thermocouple 
was  spot-welded  onto  the  test  specimen  and  protected 
by  No.  33  Sauereisen  Sealing  Cement.   Each  was 
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Figure  23:   Constraining  Equipment 
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located  0.536  inch  (1  second  of  arc  travel)  ahead 
of  the  strain  gage  mid-axis. 

c.  Instrumentation.   Strain  gages  were  connected 
into  a  Potentiometric  Circuit  (Half-Wheatstone 
Bridge) ,  balanced  and  calibrated  as  indicated 
schematically  in  Figure  24.   Thermocouples  were 
referenced  to  a  3  2°F  ice-bath  and  calibrated  as 
indicated  in  Figure  25.   Both  circuits  were  fed  into 
a  Honeywell  continuous-recording,  12-channel 
Visicorder.   When  the  raw  data  was  actually  read 
off  the  recorder  tape,  some  traces  were  delayed  or 
advanced  with  respect  to  others.   This  was  done 

to  correct  for  the  finite  difference  in  position 
along  the  weld  line  of  the  thermocouples  (1  second) 
and  strain  rosette  elements  (+  1/3  second) . 

This  produces  the  effect  of  a  simultaneous 
strain  and  temperature  reading  at  the  strain  gage 
location.   Timing  was  accomplished  by  means  of  an 
electric  stop  watch  as  well  as  a  timer  integral 
to  the  Visicorder.   Static  strain  measurements 
were  taken  with  a  Budd  Instrument  Division  Model 
P-350  strain  indicator. 

d.  Extensiometer .   Figure  26  shows  a  sketch  of 

the  extensiometer  utilized  in  Tests  5  and  6.   Strain 
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Figure  24:   Strain  Gage  Instrumentation  Circuit  (7.) 
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Figure  25:  Thermocouple  Instrumentation  Circuit  (7). 
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Figure  26:  Extensloaeter .  The  Adjusting 
Screws  are  Set  So  That  the  Spring  Holds 
Cxtensloaeter  Firaly  In  Place. 


Figure  27:   Extensloaeter  Calibration  Curve.   0.000  Inches  Represents 
an  Extensloaeter  Gap  of  3.00  Inches. 
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gages  were  located  on  both  sides  of  a  spring  steel 

ribbon  and  connected  into  adjacent  legs  of  the 

resistance  network  circuit  described  in  paragraph 

c.  above.   The  net  effect  of  such  a  set-up  is  to 

double  the  bending  sensitivity  and  null  the  stretching 

sensitivity.   Figure  27  provides  the  calibration 

curve  for  the  Visacorder  referenced  to  actual 

dimensional  changes  in  inches.   The  extensiometer  and 

its  mounting  tabs  may  be  seen  in  the  photographs 

in  Figures  28  through  31.   It  was  pulled  back  away 

from  the  torch  to  prevent  damage  when  the  arc  past 

the  mounting  tabs.   Once  the  torch  was  safely  clear, 

it  was  set  back  in  place  to  observe  weld  shrinkage. 

It  is  noted  that  the  mounting  tabs  were  cemented  in 

place  with  EPY  600  cement  vice  bolting  or  welding 

to  minimize  local  straining  of  the  base  plate. 

3.   Welding  Equipment  and  Conditions;   Welding  conditions 
are  summarized  in  Tables  7  and  8.   The  welding  machine 
utilized  was  manufactured  by  the  Linde  Division  of  Union 
Carbide  Corporation  and  consisted  of  a  type  HW-16  GMA 
torch,  a  SVI-300  power  supply,  and  associated  governor, 
carriage  and  wire  feed  mechanisms.   Travel  speed,  arc 
voltage,  amperage  and  arc  length  were  present  to  the  same 
value  for  all  experiments.   Some  fluctuations  of  amperage 
did  exist  when  the  geometry  of  the  weld  changed  from 
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bead-on-plate  to  butt.   Wire  feed  is  a  functin  of  the 
other  variables  and  was  maintained  by  the  machine  auto- 
matically*  Figures  29  through  31  are  photographs  of 
the  actual  equipment  set-up. 

E.   Experimental  Procedures: 


The  experimental  operation  is  shown  schematically  in 
Figure  32.   The  test  plate  was  instrumented  and  clamped 
into  place.   The  welding  machine  was  lined  up  with  the 
joint  and  positioned  over  the  run-off  tab  at  the  left 
end  of  the  backing  plate.   Welding  speed,  arc  voltage  and 
amperage  were  pre-set.   The  visicorder  was  actuated  and 
an  arc  was  struck  on  the  backing  plate.   As  the  welding 
torch  crossed  from  starting  tab  to  plate,  the  timer  was 
started.   The  Visicorder  output  was  marked  when  the  arc 
passed  the  center  strain  gage  location.   When  the  welding 
bead  reached  the  run-off  tab  at  the  right  end  of  the 
plate,  the  arc  was  extinguished  and  the  plate  allowed  to 
cool.   Thr  recorder  continued  to  monitor  the  gages  until 
the  plate  cooled  to  ambient  temperature.   After  cooling, 
the  plate  was  released  from  its  clamp,  and  the  new  strains 
recorded  either  with  the  recorder  or  the  Budd  instrument  ^ 


74 


GMA  Welding  Machine 


Clamp-on  Ammeter 


•  Thermocouple 


/iv  Strain  Gage  Rosette 


Figure  32:  Schematic  of  Apparatus  and  Procedure.  (7) 
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IV:       DATA   DEDUCTION 


A.  General: 

At  any  point  on  a  free  (unloaded)  surface  of  a 
solid,  it  is  necessary  to  know  three  independent  quantities 
in  order  to  specify  the  state  of  stress  completely.   These 
quantities  are  the  magnitudes  of  two  principal  stresses, 
a,,  and  a2,  and  their  directions,  (j)  or  (<J>  +  90°),  with 
respect  to  some  reference.   For  isotropic  elastic  materials, 
these  values  can  be  calculated  from  strains  measured  on 
the  surface  at  the  point  in  question,  and  since  three 
independent  quantities  are  to  be  determined,  it  will  be 
necessary  to  make  three  independent  measurements  of  strain 
(42,  pp.  G2-03) . 

B.  Strain  Calculations: 

Equation  47  describes  the  relationship  for  general 
strains  for  a  three-element  rosette  shown  in  Figure  33. 


Figure  33 
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(47) 


In  the  special  case  of  a  three-element  45°  rectangular 

rosette  shown  in  Figure  34, 

B,Y 


Figure  34 


x 


Equation  (47)  reduces  to 


111 
0  2  0 
11-1 


1 


x 


Y 


xy 


(48) 


Solving  Equation  (48)  for  e  ,  e  ,  and  Yxy/ 


1-11 
0  10 
10-1 


xy 


(49) 


The  Mohr ' s  Circle  for  this  special  case  is  shown  in  Figure  35 
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Figure  35:  Mohrfs  Circle  for  the  Rectangular  Rosette  With 
Three  Observations  of  Strain  (42,  pp.  G2-12). 
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Principal   strains   are 

ei  ■  ^-^rX)  +  \  /(ex  -  y2  +  \f      '50» 


e9   =   (  X  o  Y)   "  T  /(e   "  e  )   +   Y         (51) 
2         2       2/x    y       '  xy 


where  4>  is  defined  as  the  angle  of  one  of  the  principal 
axes  with  respect  to  the  axis  of  reference.   Mathematically, 

y 

Tan  2 <(>    =   £Y_  (52) 

e   -  e 
x    y 

In  the  literature,  there  is  considerable  confusion 
concerning  4>  and  which  principal  axis  it  defines. 
Murray  (42,  pp.  G2-11,  15)  offers  a  consistent  method 
for  dealing  with  this  problem.   If  ^is  defined  as  the 
angle  measured  (positive  in  the  anti-clockwise  direction) 
from  the  positive  OA  axis  of  the  strain  rosette  to  the 
positive  01  axis  which  corresponds  to  the  direction  of 
e, ,  the  algebraically  largest  principal  strain,  then: 

e   +  £ 

Rule  1.   When  e  >  -^-= — -  ,  then  0°  <  <t>   <  90°       (53) 

b       2  1 

£   +  £ 

Rule  2.   When  e  <   -a— = — -  ,  then  -90°  <  <j>  <  0°      (54) 

b       2 

£   +  £ 

Rule  3-   When  e,   =   -*-— — £  (55) 

b        2 

and  (a)   £   >  e   then  e  =  £   and  <t>   =  0°  (55a) 

a       C  a.  J.  -L 

or   (b)   £   <  £  ,  then  e   =  £   and  4>   =  +  90°  (55b) 
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Rules  1  through  3  are  defined  with  axis  OA  as  the 
reference.   To  convert  to  the  situation  consistent  with 
this  investigation,  axis  OA  is  +45°  (anti-clockwise)  from 
the  x-axis. 


A 


B 


Mathematical  model 


weld  line 
x   x-axis 


Geometry  utilized  in 
this  investigation 


B 
Figure  36 


*  x   =   *1   +   45° 


(56) 


The  geometry  utilized  in  this  investigation  reduces 
to  that  of  the  mathematical  model  if  the  axes  are  kept 
consistent  with  Figure  36. 

C.   Stress  Calculations: 


This  investigation  yields  mechanical  strains.   The 
automatic  temperature  compensating  feature  of  the  strain 
gage  coupled  with  the  apparent  strain  correction  removes 
the  thermal  strain  effect.   Stresses  must  be  determined 
by  total  strain. 


emech 


thermal 


-   e 


total 


(57) 
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where  T_ 


thermal     =      L      o(T)    dT  <4°> 

Jti 

For  general  two-dimensional  stress  calculations, 

a       =  5-2   (e      +  ve      )  (58) 

1  -  V    xtotal      Ytotal 


ay  "   ^~2   (£v      +   Vex      )  {59) 

y     1  -  vz    *total      xtotal 


<*i   =  ^"2   (e,      +  ve      )  (60) 

1  -  vz    xtotal      ^total 


°2   =   ^"2   (£2      +   V£l      ) 

1  -  v*    ^ total      ±total  (61) 


For  this  particular  configuration, 


r  e   -  e   =   Y  (62) 

xymech     a    c      x*total 


Txy  =   G^xy  =   2(1  +  v)   ^xy  (63) 


v  is  Poisson's  ratio.   Nominally,  v  =  .333  for  aluminum. 
v   will  probably  change  with  temperature,  but  should  not 
exceed  a  maximum  value  of  .500  in  the  weld  puddle  because 
of  volume  consideration  when  melted  (at  .500,  there  will 
be  no  volume  change) . 

Stress  calculations  are  dependent  upon  E(T)  and  a(T) 
assuming  v  is  relatively  independent  of  temperature.   These 
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values  may  be  determined  from  Table  5  in  Chapter  II. 
The  values  of  E(T)  given  assumes  that  the  value  from  a 
half -hour  soak  approximates  that  from  welding. 

D.   Determination  of  Plastic  Conditions: 

Masubuchi  (43)  theorizes  that  the  region  in  the 
vicinity  of  the  weld  has  undergone  plastic  deformation 
when  the  invarient  I  is  larger  than  the  value  of  yield 
stress,  where, 

1      =       (ox2      -     <Vy     +     °y2      +      3TXy2)V2  (64) 

This  value  was  derived  from  the  second  stress  deviator 
tensor  invariant  J~  (44,  pp.  41)  which  for  a  two-dimensional 
situation  is  reduced  to 

jr   =   1  (0    2   _  a   a        +  a  2   +   3    2}  =   1  3-2   (   } 
2     3  v  x       xy      y        xy  '     3 

The  Distortion  Energy  Theory  (Von  Mises '  Yield  Criterion) 
at  the  yield  point  in  simple  tension  (44,  pp.  75)  predicts 
that 

J2   =   |-02  (66) 

Equating  Equations  (65)  and  (66) 

or 

I  -  a  -  /"3J^  (68) 
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Therefore,  a  check  to  see  if  plastic  conditions  exist 
would  occur  is  that 
I  (T) 


a  (T) 


1.0  (69) 


It  is  important  to  note  that  any  single  component  of  the 
stress  field  in  Equation  (65)  may  be  greater  than  the 
yield  stress  in  simple  tension  without  having  yield  occur. 

Again,  the  yield  stress  is  determined  from  Table  5 
in  Chapter  II. 

E.   Data  Reduction  Program: 

Appendix  B  lists  the  basic  program  utilized  for 
rosettes  in  this  investigation.   It  is  modified  slightly 
to  suit  other  geometries  (rosettes  in  different  locations) 
or  slightly  different  welding  conditions.   Input  cards 
are  listed  as  comment  cards  in  the  main  program  and 
samples  are  shown  in  the  appendix.   In  the  case  where  all 
strains  are  zero  or  equal,  the  program  defaults  to  <f>   =90°. 
For  example,  if  at  time  =  0.0  the  strains  are  initialized 
to  zero,  the  resulting  <j>   at  time  =  0.0  is  meaningless.   In 
the  case  where  initial  strains  at  time  =  0.0  are  not  zero, 
the  program  will  reference  all  succeeding  strains  to 
their  initial  value.   Single  element  strain  gage  data 
may  be  processed  by  this  same  program  to  provide  temperature 
compensation . 
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V   RESULTS: 

A.   General  Trends: 

Figures  37  through  59  graphically  present  part  of 
the  experimental  data  tabulated  in  Appendix  C.   Table  7 
discusses  the  particular  problems  and  parameters  associated 
with  each  test.   Several  effects  are  discussed  below: 

1.  Arc  Efficiency.   Figures  37  through  40  show  the 
temperature  distributions  for  Test  No.  1  and  2,  and  also 
the  one  dimensional  computer  analysis  results  for  n  =  0.7  0 
and  0.75.   The  n  =  0.75  curve  seems  to  fit  the  data  best. 
This  efficiency  factor  includes  all  heat  loss  effects 
associated  with  the  heat  generation  and  cooling.   The 

tail  of  the  experimental  temperature  distribution  is  higher 
than  the  computer  analysis  because,  geometrically,  the 
experimental  plates  are  not  infinite  and  cool  slower  than 
the  computer  model. 

2.  Placement  of  Strain  Gages.   Figure  41  shows  the 
transverse  temperature  distribution  based  on  Figures  37 
through  40.   In  general,  the  peak  temperature  predicted 

at  transverse  locations  are  supported  by  the  testing  shown 
in  Figures  42  through  44.   Figure  41  provides  an  effective 
technique  for  positioning  strain  gages  to  prevent  damage 
or  exceeding  their  temperature  compensation  limits. 
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3.   Butt  vs.  Bead-on-plate  Characteristics .   Figures  41 

through  44  provide  temperature  distribution  data  for  the 
tests  with  strain  gages  installed  and  Figures  4  5  through 
51  provide  data  on  transverse  and  longitudinal  strain  dis- 
tributions resulting  from  these  tests.   Ideally,  a  butt 
weld  becomes  more  and  more  like  a  bead-on-plate  weld  for 
multipass  situations.   One  would  hope  for   a  two  dimensional 
model  which  could  cover  both  cases.   Unfortunately, 
Figures  43  and  44  show  a  different  initial  temperature 
distribution  for  bead-on-plate  and  butt  welds.   This  implies 
that  for  the  plate  thickness  utilized  in  this  study,  the 
two  types  of  weld  will  behave  somewhat  differently  in 
the  way  of  thermally-induced  strains,  moments,  and  distor- 
tions.  The  bead-on-plate  weld  is  not  a  two  dimensional 
situation,  but  the  butt  weld  clearly  is.   Figures  45,  46, 
48,  and  49  seem  to  agree  with  this  observation  for  strains, 
since  bottom  bead-on-plate  values  do  not  match  the  top  ones. 
The  higher  transverse  strains  for  Test  No.  4  are  noted, 
but  not  fully  explained.   The  strains  were  eratic  at  their 
peak  values  as  if  the  plate  was  slipping  in  its  constraint. 

4.   Correlation  Between  Experiments.   In  general, 
there  is  not  much  correlation  between  residual  strains  from 
experiment  to  experiment.   This  problem  has  been  noted  at 
MIT  with  previous  experiments  and  is  again  noted  here  in 
the  various  strain  distributions.   An  effort  to  find  a 


means  of  correlation  was  made  by  looking  at  the  angle 
between  the  largest  principal  strain  and  the  x-axis  (<J>  )  . 
This  data  is  shown  in  Figures  52  through  55.   There  seems 
to  be  little  correlation  in  the  residual  strain  region 
other  than  that  the  major  strains  are  either  close  to 
the  x  or  the  y  axis.   Along  the  edge  of  the  panel,  however, 
the  test  data  consistently  show  a  <j>   which  implies  loading 

A 

by  shear. 

Figures  56  through  58  provide  a  much  more  successful 
means  of  correlation.   The  ratio  of  the  invariant  described 
in  Equation  (64)  to  the  temperature  dependent  uniaxial 
yield  stress  from  Table  5  is  plotted  with  respect  to  time. 
For  residual  conditions,  there  appears  to  be  very  good 
correlation  between  the  "state  of  stress**1  composed  of 
all  stress  components  from  test  to  test.   In  addition,  when 
the  value  of  this  ratio  exceeds  1.0,  a  yielding  situation 
exists,  and  plastic  deformation  occurs. 

The  static  residual  strain  measurements  made  during 
this  study  offered  such  a  wide  and  inconsistent  variety  of 
values  that  they  were  not  included.   One  possible  reason 
for  this  is  that  the  position  of  the  strain  gages  one  inch 
away  from  the  centerline  of  the  plate  may  have  been  near 
the  region  in  Figure  1-4 ,  where  the  stress  distribution  went 
through  a  zero  point  with  a  large  slope  .   Hence,  a  minor 
variation  in  stress  could  imply  large  positive  or  negative 
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variations  in  strain. 

5.   Extensiometer .   The  extensiometer  results  are 
plotted  in  Figure  59.   Initial  measurements  vary  from 
experiment  to  experiment  because  of  the  difficulty  of 
placing  the  holding  tabs  in  the  exact  same  position. 
However,  it  is  interesting  to  note  that  for  the  two 
situations  presented  here,  once  the  weld  was  filled  in 
the  vicinity  of  the  extensiometer,  the  extensions  became 
a  measurement  of  transverse  shrinkage  and  were  similar  in 
both  tests.   Test  No.  1  utilized  two  tack  welds  at  the 
ends  of  the  plate,  and  the  arc  became  unstable  as  a  result 
of  metal  movement  in  the  center.   All  subsequent  tests 
utilized  four  equally  spaced  tacks,  so  the  plates  were 
essentially  rigid  except  in  the  immediate  region  of  the 
arc. 

B.   Accuracy  of  the  Experimental  Model  and  Instrumentation: 

1.   The  Physical  Model .   As  mentioned  previously, 
the  plate  geometry,  level  of  constraint,  and  welding  pro- 
cedures are  considered  to  be  representative  of  ship 
fabrication  processes.   However,  the  means  of  constraining 
the  panels  is  unsatisfactory  because  of  evidence  of  slipping 
Considerable  difficulty  was  experienced  in  cementing  the 
strain  gages  to  the  model.   The  best  cement  still  remains 
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EPY-600  because  it  offers  a  lower  curing  temperature. 
However,  pressure  must  be  applied  to  the  gage  while  the 
cement  is  curing,  making  this  an  extremely  critical  stage 
in  model  production. 

2 .  Data  Reduction  Calculations.   The  temperature 
compensation  has  been  previously  discussed.   However,  in 
the  strain  distributions  presented,  there  is  a  strange 
fluctuation  in  the  results  in  the  vicinity  of  the  arc 
immediately  before  the  thermal  strains  take  off.   It  is 
of  the  opinion  of  this  author  that  this  is  in  part,  if 
not  all,  due  to  temperature  compensation  techniques  and 
represents  only  an  "apparent  strain".   In  any  case,  it  is 
a  minor  effect  and  relatively  unimportant. 

3.  Instrumentation .   The  Visacorder  output  was 
scaled  initially  in  Test  No.  4  to  1200  microstrain/inch 
and  for  later  tests  at  500  to  600  microstrain/inch.   The 
reading  accuracy  is  well  within  the  5%  instrument  accuracy 
estimated  by  technicians  aiding  in  this  study. 

4.  Welding  Machine.   The  welding  machine  utilized 
in  this  study  leaves  much  to  be  desired.   The  speed  drive 
appears  not  to  give  a  uniform  arc  velocity  although  the 
time  to  span  the  full  plate  length  seems  to  be  consistent. 
The  voltmeter  and  ammeter  on  the  power  supply  have  never 
been  calibrated  and  a  clamp-on  ammeter  and  separate 


90 
voltmeter  had  to  be  utilized  as  a  result.   The  meters  on 

the  actual  welding  machine  controls  are  either  inoperative 

or  way  out  of  calibration.   There  is  no  screen  over  the 

tip  of  the  torch  which  prevents  back  splatter  from  fouling 

the  automatic  filler  metal  feed  system.   Luckily,  the  wire 

jams  which  occurred  as  a  result  were  on  the  starting  tabs 

and  the  expensive  test  plates  were  not  ruined. 
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VI   RECOMMENDATIONS : 


1.  The  small  number  of  tests  performed  in  this  study- 
should  be  expanded  into  more  tests,  ultimately  with  plates 
of  varying  thicknesses.   Data  reduction  from  the  Visa- 
corder  strip  charts  is  difficult  and  time-consuming.   If 

a  large  number  of  test  are  run,  it  becomes  almost  imperative 
that  some  form  of  direct  digital  output  be  utilized  to 
provide  punched  tapes  or  cards  compatible  with  MIT ' s  IBM 
computer  system.   Multipass  tests  likewise  require  this 
capability.   Welding  is  a  statistical  process,  so  the 
larger  the  basis  of  data  the  better  the  results  will  be. 
Obviously,  there  are  serious  questions  of  financial  and 
manpower  resources  available,  so  the  trade-off  of  more 
experiments  vs.  better  data  must  be  seriously  considered. 

2.  More  bead-on-plate  tests  which  record  top  and  bottom 
strain  and  temperature  distributions  are  needed  to  verify 
a  two  dimensional  approximation. 

3.  The  welding  machine  utilized  by  this  study  should  be 
carefully  inspected  and  the  discrepancies  noted  in  Chapter 
V  remedied,  or  another  machine  must  be  found. 

4.  A  "V"  or  "J"  weld  preparation  should  be  considered 
for  butt  welding  thicker  plates. 
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5.  Smaller  strain  gages  should  be  utilized  to  allow 

easier  application  and  closer  positioning  to  the  weld. 
This  sacrifices  the  larger  area  for  averaging  measurements, 
however.   Multiple-element  gages  should  be  of  the  stacked 
variety  to  eliminate  time  corrections  for  individual 
elements,  and  make  data  reduction  simpler.   Gages  such 
as  BLH  Electronics1  FABR-12-12SX13ET  (stacked,  three-element 
45°,  0.12  inch  gage  length),  FABX-12-12S13ET  (stacked,  two- 
element,  90°,  0.12  inch  gage  length),  or  FAER-06RB-S13ET 
(three-element,  45°,  .0625  inch  gage  length)  or  equivalent 
Micro-Measurements,  Inc.  gages  should  be  utilized. 

6.  If  time  permits,  an  independent  check  of  apparent 
strain  corrections  is  recommended  to  verify  factory  data. 
This  may  be  accomplished  by  cementing  a  gage  to  a  sample 
test  plate  and  heating  it  in  an  oven. 

7.  High  temperature  strain  gages  should  again  be  con- 
sidered to  see  if  they  are  any  more  reliable  than  reported 
previously  by  Klein  (7) . 

8.  The  test  plate  mounting  assembly  should  be  modified 
to  allow  bolting  vice  clamping  of  the  plate  boundaries. 
This  may  be  done  simply  by  drilling  and  tapping  the  bed 
plate  currently  in  use,  and  will  eliminate  any  slipping 
from  taking  place. 
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9.  Further  use  of  the  extensiometer  is  not  recommended. 
It  costs  one  channel  on  the  Visacorder  and  causes  the 
loss  of  a  rosette  (12  channels  are  available:   3  rosettes 
and  3  thermocouples,  or  2  rosettes,  2  single-element  gages, 
the  extensiometer  and  3  thermocouples) .   This  author  feels 
that  there  is  more  to  be  learned  at  present  from  the  "state 
of  stress"  than  from  shrinkage  measurements.   However,  if 
continued  use  of  the  extensiometer  is  desired,  two  more 
gages  should  be  installed  to  double  its  sensitivity. 

10.  Continued  basic  research  is  needed  to  provide  better 
physical  and  mechanical  properties  for  a  larger  variety  of 
structural  materials. 
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